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Abstract 

For a three-body system, a quantum wave function with definite £ and 
m quantum numbers may be expressed in terms of an internal wave function 
xl. which is a function of three internal coordinates. This article provides 
necessary and sufficient constraints on xi ensure that the external wave 
function is analytic. These constraints effectively amount to boundary 
conditions on xi ^^'^ derivatives at the boundary of the internal space. 
Such conditions find similarities in the (planar) two-body problem where the 
wave function (to lowest order) has the form rl"*! at the origin. We expect the 
boundary conditions to prove useful for constructing singularity free three- 
body basis sets for the case of nonvanishing angular momentum. 
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I. INTRODUCTION 



Consider a wave function for a system of three bodies that is an eigenfunction of 
and Lsz with quantum numbers I and m, respectively, where is the space-fixed orbital 
angular momentum. (The s subscript stands for "space- fixed" . ) We regard '^1^-, the "external 
wave function," as a function of the two Cartesian Jacobi vectors. It is well known 
that such a wave function can be written in the form 

'^L = Y.'^tkxl (1.1) 

k 

where the Wigner rotation matrix "D^j^. is a function of the three Euler angles and where Xki 
the "internal wave function," is a function of three internal or shape coordinates. As usual, 
k is regarded as the quantum number of the body- fixed L^. The wave function \Ef^ need 
not be an energy eigenfunction; for example, it could be an element of a basis set in terms 
of which an unknown energy eigenfunction is to be expanded. The basis sets we have in 
mind include standard orthonormal bases, hyperspherical harmonics [^0, discrete variable 
representation (DVR) bases ||6|-[TT|, and wave packet or wavelet [0 bases. 

This paper concerns boundary conditions which the internal wave function Xk iiiust 
satisfy, given that the external wave function \E'^ is a smooth function of the Cartesian 
Jacobi vectors. The boundary in question is the boundary of shape space, which consists of 
the coUinear configurations. The applications we have in mind are mainly molecular (either 
bound states of triatomic molecules or triatomic scattering problems), but the considerations 
we raise also apply to atoms or other systems of bodies with Coulomb interactions (with 
certain qualifications discussed below). We ignore spin in this paper. We consider only three- 
body systems in this paper, but an important reason for studying boundary conditions in 
three-body systems is that it is good practice for the analogous problem for systems of four 
or more bodies, which is generally more difficult and much less well understood. 

There are at least two practical reasons for being interested in boundary conditions. First, 
if one attempts to expand some unknown function in terms of a given basis, and if there are 
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boundary conditions satisfied by the unknown function wliicli are not satisfied by the basis 
functions, then the convergence will be slow. In important cases, the coefficient of the n-th 
term in the expansion will fall off either exponentially or algebraically with n, depending on 
whether the basis does or does not satisfy the required boundary conditions, respectively. 
For example, it is a bad idea to use the ordinary Legendre polynomials to expand a function 
whose 6 dependence has the boundary conditions of one of the associated Legendre functions 
for m > 0. The importance of properly treating such boundary conditions in three-body 
problems has been discussed previously by Kendrick et. al. 



Sometimes basis functions are created on the internal space simply by writing down 
some internal wave functions xi that are considered convenient, for example, distributed 
Gaussians or wave packets. From the given internal wave functions, one can construct the 
corresponding external wave functions according to Eq. (|1 . 1| ) . The question then arises, will 
these external wave functions have the same smoothness and analyticity properties as some 
unknown wave function (usually an energy eigenfunction) which one wishes to find? If not, 
the convergence will be poor. For another example, it is common practice to create internal 
basis functions by writing down the exact internal Hamiltonian, and then carving out some 
piece of it which has eigenfunctions which can be determined analytically. Again there is 
a question as to whether the basis functions created in this manner have the boundary 
conditions required of the desired unknown eigenfunction. The answer to this question 
depends in part on whether the operator created by carving out a piece of the Hamiltonian 
is itself well behaved. The analysis of this paper will show how to answer these questions. 

A second reason for being interested in boundary conditions is that numerical methods 
for solving partial differential equations on a grid must generally take careful account of 
boundary conditions, in order to guarantee reasonable accuracy and convergence. Grid and 
basis set methods are related, since grid methods implicitly involve a set of localized basis 
functions associated with the grid points. For example, DVR methods involve a basis set 
consisting of localized wave functions, resembling diffraction patterns from a narrow slit. 
The often cited "unexpected accuracy" (that is, rapid convergence) of DVR methods is 



closely related to satisfying the right boundary conditions. For example, the trapezoidal 
rule converges exponentially (rapidly) when an analytic, periodic function is integrated over 
a period, but has only power law (slow) convergence if the function is analytic but not 
periodic, or if it is integrated only over a partial period. Similarly, if the wrong DVR basis 
set is used for a given problem, the (by now expected) unexpected accuracy will be lost. 
We will have more to say about boundary conditions and rates of convergence in future 
publications, but this paper will concentrate on the boundary conditions themselves. 

We will momentarily present the principal results of our analysis of boundary conditions 
in three-body systems, but first it is well to recall some facts about two-body systems with 
rotational invariance. Thinking of energy eigenfunctions, we will speak first of the problem 
of central force motion. In three spatial dimensions, the energy eigenfunction can be written 
\E'^(r) = y^(6', 0)x£(r), where x^(r) is the internal or radial wave function, defined on the 
radial half-line < r < oo, which is the internal space. According to the standard textbook 
analysis [0, the radial wave function behaves as near r = 0. This behavior holds when 
the true potential V{r) is analytic at r = 0, but also in other cases such as that of the 
singular Coulomb potential. 

The standard analysis that produces these results proceeds by expanding the radial 
wave function in a Taylor series about r = and balancing terms on the two sides of the 
Schrodinger equation. Unfortunately, this analysis leaves the impression that the behavior 
Xe ~ of the radial wave function near r = applies only to energy eigenfunctions. 
Actually this behavior is much more general. Consider any wave function ^^^(r) which is 
an eigenfunction of and and which is analytic at r = when expressed in terms of the 
Cartesian coordinates {x,y,z). This would apply to the eigenfunctions of any rotationally 
invariant operator which is well behaved at r = 0, including Hamiltonians with central 
potentials V{r) which are analytic at r = 0. With standard assumptions about phase 
conventions the wave function can be written ^^^(r) = Y^{9, (j))xe{r), where Xi{^) is the 
radial wave function. We note that this form follows from the standard theory of rotations 
and the fact that is an eigenfunction of and Lsz', we do not invoke separation of 
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variables, since we are not necessarily separating any wave equation. Then it turns out that 
Xe{f) is analytic at r = 0, and that its Taylor series begins with the term and thereafter 
contains only the powers r^"*"^", n = 1,2, . . ., that is, every other integer power of r. Energy 
eigenfunctions in the Coulomb problem do not fit this pattern, since "^^ir) has a cusp at 
r = and is not analytic there. This is because V{r) is not analytic at r = 0. Although 
Coulomb radial wave functions xi do go as near r = 0, the Taylor series of Xi contains 
every subsequent power of r, not every other one. 

It is also worthwhile mentioning the case of two bodies in two spatial dimensions, since 
planar two-body boundary conditions bear a strong analogy to the boundary conditions 
in three-body systems. For the planar problem, an eigenfunction of Lg^ can be written 
= e*™'^Xm('")) where r = {x,y), is the azimuthal angle, and m = 0, ±1, ±2, ... is 
the quantum number of Lgz- If is analytic in the two Cartesian coordinates {x,y) at 
r = 0, then the radial wave function Xmi'f') is analytic at r = 0, its Taylor series begins with 
the r''"! term, and subsequently contains only every other power of r, rl™!"*"^", n = 1,2, . . .. 

At this point a reader who works with three-body molecular problems may wonder what 
the relevance is of two-body boundary conditions at r = 0, since r = is the two-body colli- 
sion and the collision of atoms in molecular problems does not happen at ordinary energies. 
The answer is that in three-body systems the coUinear configurations play somewhat the 
same role as the two-body collision in the planar two-body problem, insofar as boundary 
conditions are concerned. This is why it is important to know about two-body boundary 
conditions at r = 0, even for molecular problems. Collinear configurations are not necessarily 
suppressed in three-body molecular problems, and are often important. 

On the other hand, collisional configurations are important in Coulomb problems, where 
the wave function has cusp-like singularities p!6|JT7[] , again because of the nonanalyticity of 
the potential. Since this paper studies the boundary conditions on the internal wave function 
which result from the analyticity of the external wave function, and since the external wave 
function in Coulomb problems is not analytic at collisions, the analysis of this paper does not 
apply to collisional configurations in Coulomb problems. But of course Coulomb problems 
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also have collinear configurations, and our analysis does apply to these, as long as they are 
not also coUisional. We can summarize by saying that the analysis of this paper applies to 
all the important boundary conditions in three-body molecular problems, and to some of 
them (the collinear, noncollisional configurations) in Coulomb problems. 

We now return to the three-body problem, and summarize the main results of this paper, 
which relate the analyticity of the external wave function to the behavior of the internal 
wave function xi ^^e boundary of the internal space. The boundary of collinear shapes 
is a plane specified by 9 = in Smith's flSl hyperspherical coordinates, or W3 = in the 
coordinates to be introduced below. One point of this plane is the three-body collision, which 
is excluded from our analysis. At other points of this plane, we have established necessary 
and sufficient conditions on xi such that the external wave function should be analytic 
functions of the six Cartesian components of the two Jacobi vectors. If is analytic at 
a collinear configuration (not the three-body collision), then the internal wave function xi 
satisfies the following properties. First, xi is itself analytic at the collinear shape, when 
expressed in terms of certain internal coordinates to be described below. Suffice it for now 
to say that one of these coordinates, call it w^, measures the mass- weighted distance from 
the plane W73 = of collinear shapes, while the other two coordinates, Wi and W2, indicate 
where we are on this plane. 

The second property involves a modified version of the internal wave function (in 
contrast to xi)- The distinction is that k is the eigenvalue of body-fixed L^, whereas fi is 
the eigenvalue of body-fixed n ■ L, where n is the body-referred unit vector specifying the 
axis of a collinear shape. The vector n is defined on the boundary plane of the internal 
space (excluding the triple collision), and is a function of where we are on that plane. Both 
k and /i range from —i to +i, and the two internal wave functions are related by a rotation 
which maps the body z-axis into the n-axis. Then, as we shall show, it turns out that if 
is expanded in powers of W3, corresponding to movement in the internal space away from 
the boundary in the direction of increasing ws, then the first nonvanishing power is Wg^', 
and subsequently only every other power occurs in the Taylor series, ^3^'^^", n = 1, 2, . . .. 
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The converse is also true: if xi satisfies these two properties at a colhnear shape (not the 
triple collision), then the external wave function is analytic. Obviously, these boundary 
conditions are like those of the planar two-body problem at r = 0, with W3 playing the role 
of r. 

We have also proved the (plausible) fact that at configurations which are not at the 
boundary of the internal space (noncollinear configurations) , the external wave function 
is analytic if and only if xi is analytic. The stated conditions, both on and off the boundary 
plane, are independent of the convention for Euler angles or the convention for body frame, 



assuming we avoid certain body frame singularities [|19]-|2T|. 

We exclude the triple collision because there is an inevitable "string" singularity ||2^ 



of the body frame in the neighborhood of this configuration that complicates the analysis. 
This configuration is not important in three-body molecular problems, but is so for Coulomb 
problems (where the singularity of the potential is a further complication). 

The outline of this paper is as follows. Section ^ contains the main result of the paper. 
Theorem |I|, which concerns the boundary conditions satisfied by the internal wave function. 
We have stated this theorem in as nontechnical language as possible. Section |T| also contains 
a description of how the boundary conditions are to be applied and discusses several explicit 
conventions for internal coordinates and body frame. Sections |T| and |I^ are more technical, 
and are devoted to proving Theorem |I]. Section |ITl| states and proves boundary conditions 
for the simpler case of the planar two-body problem. Section |^uses the results of Sect. |Il 
to prove the conditions for the three-body problem. Section contains the conclusions. 
An Appendix collecting several facts about the representations of SO (2) is included for 
reference. 



II. THE THREE-BODY BOUNDARY CONDITIONS 
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A. Review of three-body formalism 



Before stating the boundary conditions on three-body wave functions, we review some 
necessary facts about the three-body problem. We closely follow the notation and spirit 
of Littlejohn and Reinsch 0. A three-body configuration in the center of mass frame 
is described by two (mass- weighted) Jacobi vectors r^^, a = 1,2. Jacobi vectors are a 
standard topic in n-body problems The s subscript indicates that r^o, is referred to 

the space, or laboratory, frame. For obvious reasons, we will call the space of Jacobi vectors 
"configuration space" . 

It is often convenient to specify a configuration by its shape and orientation. By the 
shape, we mean the configuration modulo physical rotations; it is parameterized by three 
rotationally invariant quantities (called internal, or shape, coordinates). We will denote 
shape coordinates in general by g^, yU = 1, 2, 3. A specific and particularly useful set of such 
coordinates is (^1,^2,^3) (henceforth called the "w-coordinates" ) defined by 

wi = rl^-rj2 = p'^cos2Qcos2<^, (2.1) 
W2 = 2rsi ■ = cos26 sin2$, (2.2) 
W3 = 2|r,iXrs2| = p^sin2e, (2.3) 

where we have expressed the ti'-coordinates both in terms of the Jacobi vectors and in terms 
of Smith's |T8|j23[| symmetric hyperspherical coordinates (p, 0,$). Here p = {r'^i + r^2)^'^^ 
is the hyperradius. The w-coordinates have been used by many researchers over the years, 
including Gronwall |2^, Smith [T^, Dragt [^, Iwai [^, Aquilanti et. al. [^], and others. 



The tf -coordinates have ranges —00 < Wi,W2 < 00, < W3 < 00 and are in one-to-one 
correspondence with three-body shapes. Thus shape space is the closed half-space containing 
the physical region of coordinate space W3 > 0. Sometimes it will also be convenient to 
consider the region of coordinate space W3 < 0, which we call the unphysical region. The 
boundary plane W3 = consists of all collinear shapes. The boundary conditions to be 
presented below occur along this plane. We define w = (wf + w| + 101)^^"^ to be the "radius" 
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in w1W2Ws-sp8.ce and note the identity, 

w = {wl + wl + wlY/^ = p\ (2.4) 

The origin w = is the three-body colhsion and is an especially singular point. 

The orientation of a configuration is defined relative to some convention for body frame. 
A body frame convention may be defined by specifying the functions ra{q), where Fq, repre- 
sents the body components of the Jacobi vectors and where q represents an arbitrary shape. 
(Equivalently, q stands for {qi,q2,qz))- Our rule is to attach an s-subscript to quantities 
referred to the space frame, and to omit this subscript for quantities referred to the body 
frame. The quantities Ya{q) can also be thought of as specifying a reference orientation for 
a given shape g, relative to which other orientations of the same shape are referred. In the 
reference orientation, the space frame is identical to the body frame. The orientation of 
the configuration is defined as the rotation matrix R G 5*0(3) which rotates the reference 
configuration into the actual configuration, 

r.a = Rra(g). (2.5) 

We discuss several common choices of body frame in Sect. [11 (J|, but for now we leave the 



body frame unspecified. 

For coUinear shapes, R is not uniquely determined by Eq. ( p.5|) and there is no unique 
body frame associated with a particular choice of reference orientation. Nevertheless, the 
functions rQ,(g) are normally well defined at collinear shapes, and the assignment of a refer- 
ence orientation for a collinear shape will prove to be a useful concept. We will have more 



to say in Sect. [II B| about the singular nature of the body frame at collinear shapes. 

We now turn to the three-body wave function and give a quick derivation of Eq. ( [1 . 1[ ) 
from angular momentum theory. The wave function \I' depends on the two Jacobi vectors 
Ysa- Rotations act on such wave functions by 

(7^(Q)vI/)(r,,) = ^(Q^r,«), (2.6) 
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where TZ is the rotation operator parameterized by the rotation matrix Q G 5*0(3), and 
T is the transpose. We consider a collection of 2£ + 1 wave functions \E'^, —i<m<£, 
with definite total angular momentum i and transforming under the action of 5*0(3) via the 
standard representation, 

^(Q'^r,.) = (7^(Q)0(r,„) = $^2^L(Q)*i(rs«), (2.7) 

k 

where Vl^{Q) is the Wigner rotation matrix of Q. We use the (active) conventions of 



Messiah |15l and Biedenharn and Louck |27]. We define the internal wave function by 



Xiiq) = ^i{raiq)). (2.8) 

The internal wave function is a multicomponent wave function which we call a "spinor." 
The spinor index is k, —i < k < i. Using Eqs. ( ^.5] ), ( |2.7| ), and ( |2.8| ) we obtain the final 
result 

^i(r.«) = Ki^r^m = J2'^'r:kmi{q). (2.9) 



The importance of Eq. ( p.9|) is that it completes the one-to-one correspondence between 
the external and internal wave functions. Equation ( p.8|) gives the internal wave function 
in terms of the external wave function while Eq. (|2.9|) gives the external wave function in 
terms of the internal wave function. One must realize, however, that all 2i + 1 components 
of xi must be specified to construct whereas only one is needed to construct xi- 
This is because the different ^E'^'s are not independent but related by raising and lowering 
operations, whereas the xi's are not. 



B. Statement of Boundary Conditions 

We now turn to the principal question addressed in this paper: What are necessary and 
sufficient conditions on the internal wave function xiil) which ensure that the external wave 
function \E'^(rsa) is analytic? Naively, one might expect the only condition to be that xiil) 
is itself analytic, but this is not sufficient as we now explain. 
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First, we review some fundamental issues regarding analyticity. Further clarification 
of these points may be found in any basic reference on differential geometry, for example 



Refs. p8| , p9| . To say that a function of several variables is analytic at a point means that 
the function agrees with its Taylor series in a neighborhood of that point. To say that a 
function defined on a smooth manifold is analytic at a point means that the function, when 
represented in a suitable choice of coordinates, is an analytic function of those coordinates 
at that point. The choice of coordinates is critical since a function analytic with respect 
to one set of coordinates may not be analytic with respect to another. Thus, whenever 
we say that a function defined on a smooth manifold is analytic, we must be careful to 
say with respect to what set of coordinates. Now, if two sets of coordinates are related to 
one another by an invertible analytic function (with analytic inverse) then these two sets of 
coordinates are said (in standard mathematical terminology) to be "compatible." A function 
analytic with respect to one set of coordinates is also analytic with respect to a compatible 
set of coordinates. Thus, the analyticity of a function is defined relative to an entire set of 
compatible coordinates. 

As an example, consider the plane M^, and take the standard {x,y) coordinates as the 
privileged coordinates defining analyticity. Polar coordinates (r, 6) are compatible with (x, y) 
everywhere except at the origin (and on a radial line). Thus, a function such as /(r, 9) = r 
which is an analytic function of the polar coordinates may still not be an analytic function at 
the origin of M?. One needs additional boundary conditions at the origin to guarantee that 
a function analytic in polar coordinates is truly analytic on R^. This example contains the 
core idea of why boundary conditions may be needed to guarantee analyticity of a function. 



We explore this example further in Sect. Ill 



When considering the three-body wave function we take the privileged set of co- 
ordinates defining analyticity to be the Jacobi coordinates r^afc, a = 1,2, = x,y,z. The 
reason for choosing these coordinates is that in practice the potential energy and the wave 
functions are typically analytic (except at collisions) with respect to these coordinates. So 
long as one only uses Jacobi coordinates, this fact is sufficient and the rest of this paper 
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could be skipped. However, as already pointed out, a shape and orientation description of 
configuration space is often advantageous and this naturally involves expressing in shape 
and orientation coordinates as in Eq. ( p.9[ ). An important fact is that shape and orientation 
coordinates are never compatible at coUinear shapes, because the rotation matrix R is not 
defined by Eq. ( ^^.51 ). This is much like the relation between rectangular and polar coordi- 
nates at the origin of the plane, discussed earlier, and it explains why the analyticity of xi 
alone is not sufficient to guarantee the analyticity of \E'^ at coUinear shapes. 

We will also be interested in the analyticity of functions defined on the internal space 
(such as the internal wave function xi)- We have found it most useful to take the w- 
coordinates as the privileged coordinate system with respect to which the analyticity of 
such functions is defined. All internal coordinate systems in common use are compatible 
with the w-coordinates at most locations in shape space. On the internal space there is 
the additional issue of what we mean by analyticity at the boundary W3 = 0. We will 
say that a function defined on the internal space is analytic at a particular point on the 
boundary if the function has an analytic continuation into the nonphysical region W3 < in 
the neighborhood of that point. 

The analyticity of the functions ra{q) requires special comment. These functions have 
singularities on certain curves ("strings") in wiW2U'3-space, which emanate from the three- 
body collision w = and go out to infinity. The location of these strings depends on the 
convention for body frame [|19]-|21[. In the following we wish to work with body frame 



conventions which remove the strings from the region of interest in the internal space. 

The reader may wonder about the analyticity of the Euler angles, or of functions of them. 
As it turns out, we never need to worry about such issues, because the main result of this 
paper, which is the establishment of necessary and sufficient conditions for the analyticity 
of \1'^, only involves the analytic properties of xi- The results we prove below are valid for 
arbitrary conventions for Euler angles. 

Turning away now from general notions of analyticity, we observe that at a coUinear 
shape, which is not the triple collision, there is a well-defined (up to sign) unit vector, 
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denoted by n, pointing along the body-referred axis of collinearity. The vector n depends on 
the position along the boundary plane and is undefined off of the plane. We therefore choose 
a convention for extending n off of the boundary plane. That is, we choose a function n(g), 
defined on shape space, which points along the collinear axis when evaluated at a collinear 
shape. When evaluated at a noncoUinear shape, we only require that n(g) lie in the plane 
spanned by ri(g) and r2(g). 

We now introduce a certain basis of spinors T^{q), — ^ < /U < ^, which are eigenspinors of 
the projection n(g) -L of the body-referred angular momentum operator L onto the collinear 
axis. These spinors are chosen to satisfy 

(n ■ L)r^ = /.r^, (2.10) 



V = ^{'^iMMk = ^MM'' (2-11) 



k 

where f represents the Hermitian conjugate. Here the /i subscript does not index the com- 
ponents of Tfj, but rather labels the spinor; the 2i + 1 components themselves are denoted 
by {Tfi)k, —i<k<i and are taken with respect to the normalized eigenbasis of Lz, using 
standard phase conventions. Notice that Eqs. ( p.lO|) and ( p.ll| ) determine T^lq) only up to 



an overall phase, which is another convention we have the freedom to choose. We define an 
alternative version of the internal wave function by 

xi = rlx' = Y.WkXl (2.12) 

k 

where is the column vector (xl^j •••) Xe)'^ ■ Since the spinors are orthonormal, Eq. ( [^.12[ ) 
may be inverted to give 

xi = J2(^,)kxi. (2.13) 

Given some region of interest in the internal space, there are conventions involved in 
choosing the coordinates q^, the functions ra{q), which specify the reference orientations, 
the function n(g), which extends the collinear axis away from the boundary, and the spinors 
T^{q). We will be particularly interested in a certain class of conventions which taken 
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together we call "valid" conventions. At a noncoUinear shape, a set of conventions is said 
to be valid if is compatible with the w-coordinates and the functions rQ,(g), n(g), and 
(r^)fe(g) are analytic. In particular, this means that the shape in question does not lie on 
a string singularity. At a coUinear shape, we still require the compatibility of the internal 
coordinates with the w-coordinates and the analyticity of Ya{q), and {j^)k{q). However, 
we further require that the boundary plane be given by 53 = 0. We also require that various 



functions be either even or odd. Specifically, 

qt,{-Wi) = g^K), A* = 1, 2, (2.14) 

gsl-ws) = -53(^3), (2.15) 

(nT„)(-g3) = (ATa)(g3), (2.16) 

(P±r«)(-g3) = -(P±r„)(g3), (2.17) 

n(-g3) = %3), (2.18) 

{rM-<l^) = ir^)k{q3), (2-19) 



where we have suppressed the dependence on wi, W2, qi, and q2 and where P±{q) denotes 
the projection operator onto the plane orthogonal to n(g). 

With the preceding setup, we state the main result of this paper. 

Theorem 1 Let a configuration (which is not the triple collision) have shape q, and 
assume valid conventions for the shape coordinates q^, the reference orientation v^iq), the 
vector n(g), and the spinor Ti^{q) in the neighborhood of q. 

(i) If q is noncoUinear, then is analytic at Vga if and only if x^n (equivalently xi) is 
analytic at q. 

(a) If q is collinear, then is analytic at if and only ifxli (equivalently xi) is analytic 
at q with the Taylor series 

00 

Xi(qi, ?2, qs) = E a^„(?i, ^2)4'^'+'", (2.20) 

n=0 
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where is an analytic function of {qi,q2). 

The most striking aspect of this theorem is that the wave function grows as gg'"' away 
from the colhnear shapes. We can provide the following heuristic physical interpretation 
of this rule. If a classical three-body system, under the influence of a smooth potential, 
approaches a collinear configuration, any angular momentum about the collinear axis will 
create a centrifugal barrier which prevents the system from reaching the collinear configu- 
ration. Quantum mechanically, the centrifugal barrier acts to suppress the wave function in 
the classically forbidden region. The more quanta of angular momentum about the collinear 
axis, the more the wave function is suppressed resulting in the gg'"' growth. This interpre- 
tation is dynamical in nature since it depends on the notion of a Hamiltonian. We stress, 
however, that the derivation of Theorem |l| will depend only on notions of analyticity and 
symmetry. 



C. Explicit Examples of Boundary Conditions 



For concreteness, we analyze the boundary conditions explicitly for several choices of valid 
conventions. The first example uses the w-coordinates and a body frame which coincides 
with the principal axes. The body-referred Jacobi vectors are given parametrically by 



r2{Wi,W2,W3) 



a + b 

a + b 
2V2 



Wl\l/2 W2 



Wi \ 1/2 



W2\_ 
Z + 



Z — 



a — 6 / \ 1/2 

2V2 V ab) 



1 _j_ ^2 



abJ 



\W2\ 



a = y/w + W3, 

b = \/w — W3. 



(2.21) 

(2.22) 
(2.23) 
(2.24) 



Since the collinear axis for a collinear shape is given by h{wi, W2) = z, we define the extension 
to noncoUinear shapes by 



h{wi,W2,W3) 



z. 



(2.25) 
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Equations ( |2.21|) and (|2.22| ) exhibit a discontinuity in the function ra{wi,W2,W3) along 



the i(73-axis. In fact the reference orientation approaches a continuum of different values, 
depending on the direction of approach. The Ws-axis forms the string singularity of the 
principal axis frame. This string consists of oblate symmetric tops, and the discontinuity 
there is a direct result of the ambiguity in the choice of the principal axes due to the 
degeneracy of the principal moments of inertia. Another singularity occurs in the principal 
axis frame, arising from the double-valued nature of the frame. Upon circling the W3 axis 
once, the principal axis frame does not return to its original value, but rather has rotated by 
71. Thus, to make the principal-axis frame single- valued requires introducing a branch cut. 
In Eqs. ( p.21|) and ( p.22|) , we have chosen this branch cut to be the two-dimensional surface 



W2 = 0, Wi > 0. Such string singularities and branch cuts are common to other choices of 



body frame as well and are discussed further in Ref. . 

We comment briefly on the validity of the conventions introduced here. The w- 
coordinates are trivially valid everywhere, with the coordinate being transverse to the 
collinear shapes. The function n(g) = z is constant and hence analytic everywhere. It obvi- 
ously satisfies Eq. ( p^.lSD as well. Away from the frame singularities discussed above, ra{q) 



is analytic, and it is straightforward to show that Eqs. ( |2.16| ) and ( [^.17] ) are satisfied. From 



Eq. (p.25|) we may take (r^)fe(g) = S^k, which is clearly analytic and satisfies Eq. ( p. 191) . 
Hence these conventions are valid everywhere except at the frame singularities. 

Away from the frame singularities and away from the boundary of shape space. Theorem 
m tells us that analyticity of xii^) ^ necessary and sufficient condition for analyticity of 
^m(''so)- In order to guarantee analyticity of \E'^(rsa) on the boundary of shape space (away 
from the positive Wi-axis where there is a frame singularity), xiil) must have the Taylor 
series 

XkiWl, W2, W3) = ^ akn{Wi,W2)w^3^^'^'', (2.26) 
n=0 

where of course akn{'Wi,W2) is analytic. 

In the next example, we keep the principal axis frame, but change coordinates from the 
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w-coordinates to the symmetric hyperspherical coordinates {p, 0) defined in Eqs. ( pTl] 



- ( |2.3| ). The coordinate B is transverse to the boundary of shape space which occurs at 
= 0. These hyperspherical coordinates are compatible everywhere except the w^-axis and 
the two-dimensional surface W2 = 0, wi > 0, where they experience a coordinate singularity. 
Note that the location of the coordinate singularities agrees exactly with the location of the 
singularities in the principal axis frame. We again choose n(g) = z and (T^)fe = Sf^k- It is 
again true that the conventions are valid everywhere except at the frame singularities. Little 
modification to the form of Eq. (|2.26|) is necessary except to change the coordinates which 



gives 

oo 

xi(p,$,e) = 5^a.„(p,$)el'=l+2". (2.27) 

n=0 

Of course, the coefficients afc„(p, $) in the above equation are different from those in 
Eq. (^. 

In the next example, we again use the ty-coordinates but choose a different body frame, 
called the zxz-frame in Ref. |^ or the BF^i frame in Ref. |T^. This frame places ri along the 
positive 2;-axis and r2 in the x2;-plane. Explicitly, the reference configuration ra{wi, 102,103) 
for the zxz-fTa.me is 

ri(wi, W2, W3) = ^{w + WiY/^z, (2.28) 
r2{w^, W2, W3) = (^3x + W2z) . (2.29) 



We again take n(g) = z and (T^)fc(q') = S^k, and it is easy to see that Eqs. ( |2.14]) - (p.l9|) 



are again satisfied. Equations ( p.28| ) and ( |2.29|) exhibit a string singularity on the negative 
ti?i-axis, which consists of shapes satisfying ri = 0. Intuitively, we explain the location of 
the string by the following observation: if ri = 0, then the orientation of r2 within the 
xz-plane is not fixed. The conventions are valid everywhere off of the string; there is no 
branch cut for the zxz-{ra.me as there was for the principal axis frame. The Taylor series 
given in Eq. ( [^.26D is again applicable for the present conventions. Of course the coefficients 



<^kn{u!i,u)2) are different here and the domain of validity is also different. 
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In the next example, we continue to use the zxz-fiame, n{q) = z, and {Tfj)k{q) = Sfj,k, 
but use a different set of hyperspherical coordinates (p, C, 0) defined by 

Wi=p^sin2C, (2.30) 
= cos 2^ cos (2-31) 
W3 = p'^ cos2Csm9. (2.32) 

These are the asymmetric hyperspherical coordinates of Smith 113123] . They are compatible 



everywhere in the physical region except on the i^i-axis, where there is a coordinate singu- 
larity. Notice that this coordinate singularity includes the singularities in the zxz-fTame, 
which occur on the negative t^i-axis. The coordinate 6 is transverse to the collinear shapes, 
which occur at both 6 = and 9 = n. We concentrate on the boundary 6 = first. Since 
the conventions are valid over the entire half- plane ^ = 0, the Taylor series 

oo 

xi{pX,0) = J2<^kn{p,C)e^'\^''' (2.33) 

n=0 

is sufficient to guarantee analyticity of the external wave function. With regards to the 
half-plane 6' = vr, by our definition the conventions are not valid there because we require 
the boundary be given by ^3 = 0. Nevertheless, by using new coordinates (p, (,6' = n — 6), 
the conventions are valid on this half-plane. If the external wave function is to be analytic 
on this half-plane, xi must satisfy a Taylor series in 6' identical in form to Eq. ( ^.33[ ) 
for 6. One set of functions which do satisfy the appropriate conditions at both 6 = 
and 6 = TV are the associated Legendre polynomials P^id). An internal wave function 
Xfc(p, C) ^) = bk{p, C)Pl{6), with 6fc(p, C) analytic, therefore lifts to an analytic external wave 
function (in the region of validity). Internal wave functions of this form arise naturally when 
constructing hyperspherical harmonics. (See, for example, Aquilanti et al. I^,^.) 

III. THE PLANAR TWO-BODY PROBLEM 



Before proving the theorem on three-body boundary conditions, we discuss the two- 
body problem in the plane. Not only is this simpler case useful practice for the three-body 
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problem, but in fact our proof of the three-body results relies on the the two-body results 
presented here. 

First, we must adapt the basic concepts and notation introduced for the three-body 
problem for use with the two-body problem. The configuration space of a two-body system, 
in the center of mass frame, is just the two-dimensional plane, the relative position of one 
body with respect to the other being denoted here as e M^. The shape of the two-body 
system depends only on the separation distance r between the bodies, and we denote the 
shape coordinate by g(r). We assign to each shape q a reference orientation r(g). The 
reference orientation r(g) is simply a point on the circle of radius r{q) centered at the origin 
of M^. This fact makes the reference orientations much easier to visualize here than for the 
three-body problem; the reference orientations all lie on a curve beginning at the origin and 
intersecting each concentric circle once as it moves out to infinity. An arbitrary configuration 
Tg is given by 

r, = Rr(g), (3.1) 

where R e SO {2) denotes the orientation of the system. We will denote the rotation angle 
of R by e. 

As an example, one could choose the internal coordinate to be q{r) = r and the reference 
orientation to be 

r(r) = rx. (3.2) 

This choice produces shape and orientation coordinates (r, 6) which are the usual polar 
coordinates on M?. 

We now discuss the two-body wave function ^. The action of Q e 5'0(2) on is 

(7^(Q)*)(r,) = *(Q^r,). (3.3) 

We denote by a wave function which transforms according to the irrep m of 5'0(2), 

^miQ^Vs) = (7e(Q)*„)(r,) = e-'^'^^mivs), (3.4) 
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where a is the rotation angle of Q. The internal wave function Xmifl) is defined by 

= ^m(r(g))- (3.5) 



Given Xmifi) we may recover the external wave function \&m(rs) with the aid of Eqs. ( pTTI 
O, and 



^m(r.) = ^„^(Rr(g)) = e^\m{q). (3.6) 



Equations (|3.3|) - ( p.6|) are obviously analogous to Eqs. (^^ ) - ( p.9|) for the three-body 
problem. 

We take the Cartesian coordinates (r^^:, "'^sj;) (that is, the usual space-fixed {x,y)) as 
the privileged coordinates for defining analyticity of functions on configuration space. We 
take the radial coordinate r as the privileged coordinate for defining analyticity of functions 
on shape space. Since r is a positive quantity, we must give special consideration to the 
boundary r = 0. A function on shape space is said to be analytic at r = if it can be 
analytically continued (as a function of r) into the region r < 0. 

For the three-body problem, there were four different conventions that had to be specified. 
For the two-body problem, we need only specify two conventions: the shape coordinate q and 
the reference orientation r(g); there is no analog of n(g) or T^(q') for the two-body problem. 
Away from the two-body collision, these two conventions are said to be "valid" if the shape 
coordinate q is compatible with r and if r{q) is analytic. At the two-body collision, q must 
still be compatible with r and r{q) must still be analytic. However, we also require that 
q = coincide with the two-body collision and that the following conditions be met 

q(-r) = -q{r), (3.7) 
r(-g) = -r(g). (3.8) 

We now state the following two-body theorem, which is analogous to Theorem ^ for the 
three-body problem. 
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Theorem 2 Let a configuration have shape q, and assume valid conventions for the shape 
coordinate q and the reference orientation r(g) in the neighborhood of q. 

(i) If q ^ 0, then is analytic at if and only if Xm is analytic at q. 

(a) If q = 0, then \E'm is analytic at Tg = if and only if Xm is analytic at with the Taylor 

series 



where the amn o^re constant complex coefficients. 
Proof 

For the entirety of this proof, when we say that a function of either the Cartesian 
coordinates or shape coordinate is analytic, we mean only that it is locally analytic at 
the specific points, or q respectively, mentioned in the statement of the theorem. 

Equation (|3.4|) shows that if \E'm is analytic at an arbitrary r^, then \E'm is analytic at any 
other orientation Q^^r^ with Q G 5*0(2) arbitrary. We therefore assume without loss that 
the specific configuration in the statement of the Theorem is the reference orientation r. 

(i) Assume g 7^ 0. Assume \l'm(rsa) is analytic. From Eq. ( |3.5| ), the fact that r(g) 
is analytic, and the fact that the composition of analytic functions is analytic, Xm{q) is 
analytic. 

Next assume Xm{(l) is analytic. By the assumption of compatibility, q{r) is analytic. 
Furthermore, r(vs) = (r^^ + r'^yY^'^ is an analytic function of the Cartesian coordinates. 
Hence, q{rs) = q{r{rs)) is analytic. Turning to the rotation matrix R in Eq. (|3.1| ), its 
rotation angle 6 is given by 6{rs,r) = arcsin (z ■ (rXr^/r^)), which is an analytic function 
of and r. (We only consider 6 in the range — 7r/2 < 9 < -k /2 since we are only checking 
for analyticity at the reference orientation 9 = 0.) Furthermore, r{rs) = r(g(rs)) is analytic 
since r(g) is analytic by the validity assumption and q{rs) was shown to be analytic above. 
Hence, 



00 




(3.9) 



n=0 



^(r,) = 0(r„r(r,)) 



(3.10) 
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is analytic. Since exp{im6) is an analytic function of 6, we conclude that 

vl/^(r,) = e'"''^''^hM^s)) (3.11) 

is an analytic function of the Cartesian coordinates. 

(ii) Assume q = 0. Assume is analytic. By the same argument as in case (i), 

Xmig) is analytic at g = 0. To prove Eq. ( p. 91) , we differentiate Eq. ( p.4|) d times with respect 
to r^, 

E Qn^^-Q^.^^ ■ ■ ■ g^^m) (Q^r.) = e— (-^ ■ ■ ■ (r.), (3.12) 

u, ski ^' ska / sji ^' sja / 



ki...k 



where Qjk, j, k = x,y, are the components of Q. Evaluating Eq. (|3.12|) at = produces 



Equation ( p.l3| ) shows that the rank d tensor {d/drsji...d/drsja'^m){0) transforms under the 
completely symmetrized action of S0{2) as an irrep of 50(2) labeled by m. (See Appendix.) 
The decomposition of the fully symmetrized action of 5*0(2) on rank d tensors decomposes 
into irreps as shown in Eq. ([Al|) . If m does not label one of the irreps included in Eq. ( [AID , 
that is, if m d,d — 2, —{d — 2),—d, then {d/drsji...d/drsjj^'^m){0) = 0. Consequently, 
the chain rule and Eq. (p.5|) show that 



when m ^ d,d — 2, —{d — 2), — d. This shows that the appropriate Taylor coefficients 
vanish in order to produce the Taylor series in Eq. (p.9|) . 

Now assume that Xm is analytic at g = with Taylor series shown in Eq. ( p.9|) . We 
first complete the proof of Theorem ^ assuming the internal coordinate q{r) = r and the 
reference configuration r(r) = rx. We explicitly construct a function \E'm(rs) by 

oo 

^m{rs) = Y. E W^'") ...... (3-15) 



n = A.l...fc|„|+2n 
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where the amn are the same Taylor coefficients as in Eq. ( |3.9|) and the tl^'"''^"' are the rank 
\m\ + 2n tensors defined by Eq. (|A2|) . Clearly the transformation property Eq. (|A3|) of the 
^|m|+2n g]^Q^g ^^g^^ satisfies Eq. ( p.4|) . Hence, we may apply the same analysis to as 
we have for "^m- In particular, is uniquely determined via Eq. ( |3.6| ) by the internal wave 
function Xm defined by Eq. ( |3.5| ) 

oo 

^™(rx) = ^a^„rH+2n, (3.I6) 

n=0 



Xm[r) 



where we have used Eqs. (|3.15|) and (|A4|) . Since Xm and Xm have the same Taylor series. 



Xm = Xm- Furthermore, the unique correspondence between internal and external wave 



functions guarantees that '^m = ^m- From Eq. (|3.15|) , we see that = is analytic at 
by construction. This completes the proof of Theorem § for the specific conventions chosen 
above. 

We mention two noteworthy special cases of Theorem 0, which the above analysis has 
now proven. First, if /(r) is an even analytic function, then /(r^) = f{r{rs)) is an analytic 
function of the Cartesian coordinates. (This fact is quite trivial to prove from scratch 
by simply considering the Taylor series of the two functions.) Second, if /(r) is an odd 
analytic function, then /(fs) = /(r, 6) = exp{i6) f (r) is an analytic function of the Cartesian 
coordinates, where (r, 9) are the standard polar coordinates. 

We now assume arbitrary valid conventions r(g) and q{r) for the reference orientation and 
shape coordinate respectively. To complete the proof of Theorem |^ for these conventions, 
we first define and discuss three useful functions 'F^Tg), G{rs), and H(rs) related to these 
conventions. By the validity assumptions, r(r) is an odd analytic function in the neighbor- 
hood of r = 0, and hence r(r) = r(r)/r is an even analytic function in the neighborhood of 
r = 0. Define the matrix-valued function F(r) by 

F(r) = xr^(r) + y(zxr(r))^. (3.17) 

Notice that F(r) is the unique matrix in 5*0(2) satisfying 

F(r)r(r)=x. (3.18) 
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Since r(r) is an even analytic function, we see from Eq. ( p.l7| ) that F(r) is an even analytic 
function. Hence, from the first of the two special cases of Theorem ^ mentioned above, 
F(rs) = F(r(rs)) is also analytic. Notice of course that F(rs) is invariant under all rotations 
Q e S0{2), 

F(Qr,) = F(r,). (3.19) 

By the validity assumptions, q{r) is an odd analytic function of r. Hence, from the 
second of the two special cases of Theorem ^ mentioned above, the function 

G(r,) = e''q{r) (3.20) 

is an analytic function of the Cartesian coordinates. (Here {6, r) are the standard polar 
coordinates in M^.) We define a vector-valued version of G{rs), denoted G(rs), by 

G(r,) = Re(G'(r,))x + ^(^(r,))^. = g(r)r,. (3.21) 

Obviously, G(rs) is also an analytic function of the Cartesian coordinates. Recall that the 
compatibility of the coordinate q guarantees that q{r) has an inverse, which we denote by 
q^^, that is, q~^{q{r)) = r. The function G(r5) therefore has an inverse given by 

G-\r,) = q-\r)r,, (3.22) 



as may be verified by inserting this formula into Eq. (|3.21| ). Equations ( p.21| ) and (|3.22|) 
easily admit the following results 

G(Qr,) = QG(r,), (3.23) 
G-\Qrs) = QG-\rs), (3.24) 

where Q G SO (2) is arbitrary. 

We define a new function H(rs) by 

H(r,) = F(r,)G(r,). (3.25) 
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This function has several important properties. First, since both G(rs) and F(rs) are ana- 
lytic, H(rs) is analytic. Second, 



H(r(g)) = F(r(g))G(r(g))=gx, 



(3.26) 



which follows from Eqs. ( |3.18D and ( |3.21| ). A third fact is that H is invertible. This fact 
requires more work to prove, which we do by explicitly constructing H^^. Let r'^ = H(rs). 
Then, 



r. = G-^(F-i(r.K) = F-^(r.)G-^(r:), 



(3.27) 



where we use the definition of H and the fact that G is invertible in the first equality and 
the second equality follows from Eq. ( p.24|) . Now, the magnitudes of and are related 

by 



|F(r,)G(r,)| = |G(r, 



g(|rs|), 



(3.28) 



where we have used Eq. ( |3.21| ) in the final equality. Turning this relation around and using 
Eq. (i;^ we have. 



r.|=g-^(|r:i) = |G-^(r;)|. 



(3.29) 



As witnessed by Eq. (|3.19| ), F(rs) depends only on the magnitude of its argument, and hence 
Eqs. (13.271 ) and (|3.29|) combine to produce 



r, = F-HG-^(r'J)G-^(r'J 



(3.30) 



which gives in terms of r'^. Hence, 



H-^(r;) = F-HG-^(r;))G-^(r;) 



Another fact regarding H is that if Q G 5*0(2) is arbitrary, then 



H(Qr,) = F(Qr,)G(Qr,) = F(r,)QG(r,) = QF(r,)G(r,) = QH(r,), 



(3.31) 



(3.32) 
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where the second equahty follows from Eqs. (p.l9|) and (|3.23|) and the third equality follows 
from the commutativity of the group 5*0(2). From Eq. ( |3.32|) follows an analogous identity 
for H-i 



H-i(Qr,) = QH-i(r,). (3.33) 

We now have the proper background to complete the proof. Since H is invertible, we 
introduce the function \E'm by 

§„(r,) = ^„(H-i(r,)), (3.34) 

^„(r,) = §„(H(r,)). (3.35) 

We see from Eq. ( p.33| ) that satisfies Eq. ( p.4| ) since satisfies Eq. (|3.4| ). We define 
the internal wave function Xm{(l) using the old convention r(g) = gx, 

Xm{q) = ^m(gx). (3.36) 

Using the new convention v(q), we have Xm{(l) given by 

Xm{q) = ^„^(r(g)) = ^™(H(r(g))) = ^„(gx), (3.37) 



where the second equality follows from Eq. (|3.35|) and the third from Eq. (|3.26| ). Thus, the 



functional form of Xmiq) and Xm{<l) are identical. Since we have assumed that Xm is analytic 
at g = with the Taylor series in Eq. ( |3.9|) , Xm is also analytic with the identical Taylor 
series. Now since we have proved Theorem ^ for the old conventions used to define Xm{(l), we 
have that \E'm(i"s) is an analytic function of the Cartesian coordinates. Furthermore, since 
H(rs) is analytic, Eq. ( p.35|) implies that \l/m(rs) is analytic. QSV. 

IV. PROOF OF THREE-BODY BOUNDARY CONDITIONS 

In this section, we prove Theorem 
Proof 
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For the entirety of this proof, when we say that a function of either the Jacobi coordinates 
or shape coordinates is analytic, we mean only that it is locally analytic at the specific points, 
Tsa or q respectively, mentioned in the theorem. 

A common fact we will use several times is that the Wigner matrices ^^mA;(Q) analytic 
functions of the rotation matrices Q G SO (3). Equation (^.7]) thus shows that if \E'|„ is 



analytic at an arbitrary r^Q,, then is analytic at any other orientation Q-^r^Q, with Q G 
SO (3) arbitrary. We therefore assume without loss that the specific configuration Vga in the 
statement of the Theorem is the reference orientation r^. 

(i) Assume q is noncollinear. The proof here is a straightforward generalization of the 
proof of part (i) of Theorem ^j. First assume \l/(„(rsQ,) is analytic. From Eqs. ( p.8|) and 
( p.l2| ), the fact that both ra{q) and {T^)k{<l) are analytic, and the fact that the composition 
of analytic functions is analytic, and xHq) are both analytic. 



Next assume xfiiQ) is analytic. From Eq. ( p.l3|) and the fact that {Ty)k{q) is analytic, 
Xk{<i) is also analytic. The validity assumption guarantees that q{w^) is analytic. From 
Eqs. ( ^.11 ) - ( |2.tj| ) it is evident that the functions w^irsa) are themselves analytic. Therefore, 
qi^sa) = qiw^ivsa)) is analytic. 

We now focus on the orientation matrix R in Eq. ( ^.51 ) which may be expressed in terms 
of the vectors Vga and as 

R(r,„,r„) = ^—^ [r,iv[ + r,2V^ + (r,iXr,2)(riXr2)^] , (4.1) 

|riXr2| 

vi = -r2X(r2Xri), (4.2) 

V2 = -riX(riXr2). (4.3) 



To confirm the above expression for R, we need only verify Eq. (|2.5|) and show that 
R(riXr2) = rsiXrs2) both of which follow easily from the simple identities 

v„ -r^. = |rlXr2|^ (4.4) 

v„ • (riXr2) = 0, (4.5) 

vi ■ r2 = V2 ■ ri = 0. (4.6) 

27 



Since ri Xr2 does not vanish, it is clear that R(rsa, fa) is an analytic function of the vectors 
and Fq,. Since both ra{q) and q{rsa) are analytic, the function ra{rsa) = ra{q{rsa)) 
is analytic. Hence, R(rsQ,) = R(rsa, rQ,(rsQ)) is also analytic. Recalling that the Wigner 
matrices are analytic functions of R and that q{rsa) is analytic, we see that 

Kirs.) = 5^Pi\(R(r,„))xi(g(r,J) (4.7) 

k 

is an analytic function of the Jacobi coordinates. 

(ii) Assume q is collinear. We define a new set of coordinates on configuration space 
consisting of two shape coordinates, two orientation coordinates, and two coordinates con- 
structed from the one remaining shape coordinate and the one remaining orientation coordi- 
nate. We call these coordinates the "mixed coordinates". The majority of the remainder of 
the proof will be dedicated to defining the mixed coordinates and proving the most impor- 
tant fact about them, that they are compatible with the Jacobi coordinates at the collinear 
configuration. 

First, we discuss the parameterization of the rotation matrix R in Eq. ( |2.5| ). Consider 
an arbitrary unit vector e lying in the northern hemisphere (excluding the equator), with 
n(g) as the north pole. We define U(g, e) G SO (3) to be the unique rotation such that 

U(g,e)n(g) = e (4.8) 

and such that its rotation axis lies on the equator. Specifically, the rotation axis of U lies in 
the direction of 

a(g,e) = n(g)xe, (4.9) 

and the rotation angle of U has the value arcsin |a|. Thus, U(g, e) is exphcitly given by 

U(g,e) =exp[/(|a(g,e)|)(a(g,e)x)], (4.10) 
arcsin a; 

f{x) = , (4.11) 

X 

where we have introduced the notation ax for the 3x3 matrix which maps an arbitrary 
vector V into ax v. The purpose of Eqs. ( [4.10|) and ([4.11|) is to demonstrate the analyticity 
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of U(g,e). Observe that arcsin(x) is an odd analytic function on the interval (—1, 1), and 
hence, f{x) is an even analytic function on (—1, 1). This in turn implies that /(a) = /(|a|) 
is an analytic function of a. From Eq. ( [4.9| ) and the fact that n(g) is analytic at = 0, 
a(g, e) is analytic at gs = and at all e in the northern hemisphere defined by the north pole 
n(g3 = 0). The analyticity of the exponential function permits the following final statement. 
The function U(g, e) expressed in Eq. (|4.10|) is an analytic function at ^3 = and at all e in 
the northern hemisphere with n(g3 = 0) as the north pole. 

For an arbitrary rotation angle 6 we define \/{q,6) to be the rotation by 6 about n(g). 
Explicitly, this rotation is given by 

Viq,9) = n(g)n^(g) + [sin^(n(g)x) + cos0]Px(g). (4.12) 

We now express the rotation matrix R in Eq. (|2.5| ) by 

R(e,^) = U(g,e)V(g,^), (4.13) 

where we have parameterized R by the quantities e and 6. 

We assume without loss of generality that n(g3 = 0) = z. (This assumption simply 
amounts to a judicious choice of space frame.) We may thus define W(g) G S0{3) in the 
neighborhood of 53 = to be the unique matrix such that 

W(g)z = h{q) (4.14) 

and such that its axis of rotation lies in the xy-plane. In fact, we see from Eq. ([4.8|) that 

W(g) = U^(g,z). (4.15) 

The analyticity property of U implies that W(g) is analytic at q^ = 0. We use W(g) to define 
an orthonormal frame ni{q),i = 1, 2, 3, by 

ni(g) = W(g)x, (4.16) 
n2(g) = W(g)y, (4.17) 
n3(9) = W(g)z = n(g). (4.18) 
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The functions ni{q) are of course analytic. 

The unit vector e is determined by only two of its components, of which we choose 
ei = e ■ fii and 62 = e ■ n2. The third component 63(61,62) — v"^! ~r ^2 
function of the other two in the northern hemisphere. The function 



el + elY^'^ is an analytic 



e(g, 61, 62) = eini(g) + e2n2(g) + 63(61, 62)^3(5) (4-19) 

is of course analytic. 

We next define a pair of vectors Sa{q,9), a = 1, 2, by 

s^{q,e) = V(g,^)r„(g) = (n ■ r,)(g)n(g) + [sin^(n(g)x) + cos^] (P±r,)(g), (4.20) 



where we have used Eq. ( [4.121 ). We introduce two new variables Ui and U2 by 



Ml = qs cos 9, 
U2 = qz sin 6. 



(4.21) 
(4.22) 



We define u = q^, which is convenient notation since when u is positive, it is the radial coor- 
dinate in 'UiM2-space, that is, u = ('Ui+-U2)^''^- The vectors Sq, are conveniently parameterized 
by the new coordinates 

[P±r,)(u) 



s„(ui, U2) = (n ■ ra){u)h{u) + [u2(n(M) x) + Ui] 



u 



(4.23) 



where we have dropped the explicit dependence on qi and q2- Since f{u) is an even analytic 
function of m, f{ui,U2) = /((^i + ulY^"^) is analytic in ui and U2. This fact, together with 
Eqs. ( p.l6| ) - ( p.l8|) , shows that n, h-Ta, {P±ra)/u, and hence Sq are all analytic functions of 
Ml and U2. Reintroducing the explicit dependence on qi and q2, we see that SQ,(gi, q2, ui, U2) 
is analytic. 

Using similar arguments as above, from Eqs. ( [4.9| ) and ( [4.10| ) we see that both 
a(gi, g2, ^1, M2, e = z) and U(gi, g2, mi, ^^2, e = z) are analytic. Hence, Eq. ( [4.15| ) shows 
that W(gi, g2, Ml, M2) is analytic, from which follows, using Eqs. ( |4.16| ) - ( |4.18| ), that the 
ni(gi, g2, Ml, ^2), i = 1,2,3, are analytic. Hence, Eq. (|4.19| ) shows that e(gi, g2, mi, ^2, ei, 62) 
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is analytic. From this result we also have, using Eqs. ( |4.9| ) and ( [4.10| ), that 

a{qi,q2,Ui,U2, 61,62) = a{qi,q2,Ui,U2,e{qi,q2,Ui,U2, 61,62)) and U(gi, ^2, Mi, ^2, ei, 62) = 
U(gi, q2, ui,U2, e{qi, q2, Ui, U2, 61,62)) are both analytic. 

We define the mixed coordinates to be the variables (gi, q2,ui,U2, 61,62). We can express 
r^Q, in terms of the mixed coordinates as follows 



rsa(gi, g2, Ml, U2, 61, 62) = Rra = UVr^ = U(gi, q2, ui, U2, ei, 62)8^(^1, q2, ui, U2). (4.24) 

where we have used Eqs. ( p.5|) , (|4.13|) , and ([4.20|) . Since we have already shown that both 
U(gi, g2, Wi, M2, ei, 62) and Sa{qi,q2,Ui,U2) are analytic, we see that the Jacobi coordinates 
are analytic functions of the mixed coordinates. 

To show compatibility of the mixed coordinates with the Jacobi coordinates we need 
now only show that the mixed coordinates are analytic functions of the Jacobi coordinates. 
First, since the internal coordinates are valid coordinates, they are analytic functions 
of the w-coordinates. By inspecting Eqs. ( |2.1| ) and (|2.2|), we see that Wi{rsa) and W2(rso) 
are analytic, even at a coUinear shape. However, because of the absolute value in Eq. ( f^.3[ ) 
the same can not be said of W3(rs„). However, the function w'^{rsa) is analytic. Thus, 
we have the following lemma: Any analytic function of (^1,^2,^3) may be viewed, via 
composition, as an analytic function of r^o. Using this lemma and noting that qi{w^) and 
q2{w^) are analytic functions which by Eq. ( p.l4|) are also even in w^, we see that the two 
mixed coordinates gi(r<(„) and q2{T^sa) are both analytic functions of the Jacobi coordinates. 

Note that the coordinate qzi^i^sa) is not analytic. However, since q'^{w^) is analytic and by 
Eq. ( p.l5| ) also even in w^, g|(r<ja) is analytic. This fact allows us to generalize our previous 



lemma regarding the w-coordinates to arbitrary q. As we will have frequent need of this 
more general lemma, we record it below. 

Lemma 1 If f{q) is an analytic function which is even in qs, then f{rsa) = fiQi^^^sa)) is 
analytic. 



This lemma, together with Eqs. ( |2.16D - ( |2.1(j| ), proves the analyticity of the following 



functions: (fq, ■ n){rsa), |P±i"aP(i"sa), ''"^(I'sa), n(rsQ,). Furthermore, Eqs. (|4^ ) and ([4.1 0|) 
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show that U is even in gs and hence U(rso,e) = U(g(rsa),e) is analytic at the colhnear 
configuration and at all e in the northern hemisphere with n(g3 = 0) at the north pole. 
Equation ( |4.15| ) proves the analyticity of\N{rsa) = W(g(rsQ,)) and hence Eqs. ( [4.16| ) - ( |4.18|) 
prove the analyticity of the basis vectors hi{rsa) = ni(g(rsQ,)), i = 1,2,3. 
The vector e is determined by the equation 



e = Rn, 



(4.25) 



which follows from Eqs. ( [4.12| ), (|4.13|) and ( [4.8|) . We observed earlier that away from a 
colhnear shape, the matrix R is given by Eq. ( |4.1|) , which results in the following formula 
for e 



e = z/i(g)r,i + z/2(g)r,2, 
v„ ■ n 



(4.26) 
(4.27) 



where we have used Eq. ( [4.4|) and the fact that 



(riXr2) ■ n = 



(4.28) 



since n is assumed to lie in the plane spanned by ri and r2. We will now show that Eq. ( [4.26[ ) 
is valid not only at noncollinear configurations, but at colhnear configurations as well. More 
specifically, we will prove that e(rso) is analytic at colhnear configurations. 

Considering Eq. ( [4.26|) in the neighborhood of a colhnear configuration, we first assume 
that neither ri nor r2 vanishes. Obviously, to show that e(rsa) is analytic, we need only show 
that z/Q,(r5Q,) = t'o(q'(rsa)) is analytic. From Lemma ^ this amounts to showing that t'a('?) 
is analytic and even in q^. That i^aiQ) is even in is easily proved by inserting Eqs. ( |4.2|) 
and ( [4.3| ) into Eq. ( ^.27| ) and using Eqs. ( |2.16| ) and ( |2.17] ). To show that Uail) is analytic 
requires the following observation: The ratio of two functions which are both analytic at a 
given point is also analytic at that point provided that the ratio is not infinite. (This fact 
is easily proved by considering the Taylor series of the two functions.) As qs tends toward 
0, by our previous assumption, does not vanish, nor does n. The quantity Vq,, however, 
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does vanish, but since it appears to first order in both the numerator and denominator of 
Eq. ( ^4.27| ), the ratio is finite at ^3 = 0. Thus, 1^0,(5) is analytic, and hence //^(rsQ,) and 
consequently e(rsa) are analytic as well. 

If one of the Jacobi vectors, say Tsi, vanishes at the collinear configuration, the preceding 
analysis must be modified due to the appearance of ri in the denominator of Eq. ( [4. 27] ). In 
this case, we define a new set of vectors 

t,i = r,i + r,2, (4.29) 
ts2 = r^i - rs2, (4.30) 

and their counterparts ta in the body frame which satisfy 

tsa = Rta, (4.31) 



analogous to Eq. ( p.5| ). Notice that neither ti nor t2 vanishes at the collinear shape since 
this could only occur if r2 were also to vanish, which only occurs for the triple collision. 
Now, Eqs. ( [4.1[ ) - (|i^), ( [4.26| ), and ( [4.27| ) are all valid if one replaces and by tga 
and ta respectively. We may repeat the same line of reasoning as above to show that the 
new Va{(l)i with Fq, replaced by t^, is analytic and even in ^3. Thus from Lemma ^ the 
new ^a{^sa) and hence e(rsQ,) are both analytic. Since hi{Ysa) and n2(rsQ) are analytic, we 
finally find that the two mixed coordinates ei{Ysa) = ^{rga) ■ nj(rso), i = 1,2, are analytic 
functions of the Jacobi coordinates. 

We turn now to the final two mixed coordinates ui and U2- They are defined by 
Eqs. (|4.21| ) and ([4.22|) , but may also be expressed as 



Ml = g3 cos 9 = ill - (g3VP^)ni = fii ■ U^(g3RPx)ni = fii ■ U^Mfii, (4.32) 

U2 = qs sin e = h2- U^Mni, (4.33) 

where 

M = gsRP^. (4.34) 
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The second equality in Eq. ([4.32|) follows from Eq. ([4.12|) and the orthogonality of rii and 
n. The third equality follows from Eq. ( [4.13| ). 

Recall that U(rsa,e) and e(rsQ,) are analytic. Hence, U(rsa) = U(rsa, e(rsa)) is also 
analytic. We also recall that the ni{rsa), = 1, 2, 3, are analytic. Thus, inspecting Eqs. ( |4.32| ) 
and ( |4.33| ), we need only show that M(rsQ) is analytic. The matrix M may be written with 
the aid of Eq. as 



M(g,r, 



Pa 
P3 



rsiPi(g) +rs2P2(g) + (rsiXr^2)pr(g), 

, a = 1,2, 

Va ■ r« 

g3P±(riXr2) 

Iri Xr2P 



(4.35) 
(4.36) 

(4.37) 



We first assume that ri and r2 do not vanish at the coUinear shape. To prove that M(rsQ) = 
M(g(rsa), Tsa) is analytic, we need only show that Pa{''i^sa), a = 1, 2, and P3(rsa) are analytic. 
From Lemma |l] we must therefore show that Pa{(l), oi = 1,2, and pslq) are analytic and even 
in q^. Using the definitions Eqs. ( |4.2|) , ( [1.3D , (|4.36|) , and ( [4. 371) and the validity conditions 
Eqs. ( p.l6| ) and (|2.17|) , it is straightforward to prove that PQ-(rso), a = 1,2, and P3(rsQ) 
are even in ^3. As with our proof of the analyticity of z/Q,(g), to prove that Pa{q), a = 1,2, 
is analytic all we need to do is show that the ratio in Eq. ( |4.36|) is not infinite when q^ 
tends toward 0. (Both the numerator and denominator are analytic functions of q.) This 
fact follows readily from our assumption that does not vanish and from the fact that 
Vq, which does vanish at 53 = 0, appears linearly in both the numerator and denominator. 
Proving that P3{q) is analytic requires a bit more care. Since both the numerator and the 
denominator in Eq. ( [4.37| ) are analytic functions, we again need only verify that P3(g) does 
not blow up at g3 = 0, whence 



lim 03(^3) 



lim 

93^0 



2g3 



lim 

93^0 



ri Xr2 

|riXr2| 



(^3 



(4.38) 



where we have used the definition Eq. ( |2.3| ) of w^- We now note 

93 1 dqs 



lim 



^3-0^^3(93) 1^(0) 



dw' 



■(0), 



(4.39) 
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which cannot be infinite since gs is an analytic function of w^. (We have used I'Hospital's 
rule in the second step of Eq. ( [4.39| ).) The second limit in Eq. ( [4.38| ) cannot be infinite 



either, since P_l is well-defined at gs = and ri Xr2/|ri Xr2| is a unit vector. Thus, ps^q) is 
analytic at qs = 0. Hence, we have shown the analyticity of Pai^i^sa), « = 1,2, and P3(rsQ,) 
from which follows that analyticity of M(rsQ) and Ui{Ysa)ii = 1,2. If the assumption that 
ri and r2 do not vanish at the coUinear shape proves to be false, then we can prove the 
analyticity of M(rsQ) by applying the same trick used earlier of defining the vectors tga- We 
omit the straightforward details of completing this argument. This finishes the proof of the 
compatibility of the mixed coordinates with the Jacobi coordinates. 

We turn now to the eigenspinors satisfying Eqs. (p.lO|) and (|2.11|) . We denote by 



{Lsi)mm', i = x,y,z, m,m' = —i,...,£, the components of the space-referred angular mo- 
mentum operator Lsi with respect to the Lsz eigenbasis. The matrices Lsi transform under 
a rotation Q G SO (3) as a vector operator, 

J2 Q^Lsj = V\Q)L,,V'\Q). (4.40) 

j 

We denote by {Li)kki, i = x,y,z, k,k' = —£,...,£, the components of the body-referred 
angular momentum operator Lj with respect to the eigenbasis. The components of the 
two operators Lsi and Li, with respect to their respective bases, are related by {Li)kk' = 
{Lsi)k'k = (-^si)fcfc'- (S^^ 0; Sect. 4.H for a derivation.) Thus, the matrices Li satisfy 

J2QnLj = V'*{Q)LiV'^{Q). (4.41) 

j 

In view of Eqs. ( |4.14| ) and ( [4.41| ), we have 



n-L = z-W^L = P^*(W)L,P^'^(W). (4.42) 

From Eq. ( [4.42|) we see that can be expressed as 

{r,k = e-'-Dl, (W^) , (4.43) 

where a{q) is a phase factor which must be analytic and even in ^3 in order to guarantee 
the same properties for (r^)jt(g). This result, together with Eqs. (|2.12|) and (p.l3|) , yields 
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k 

We may use Eq. ( [4.45| ) to rewrite Eq. ( p.9|) as 



(4.44) 
(4.45) 



= e'" Yl ^r:,imxi = ^^i(UVW)x^, (4.46) 

where we have used Eq. ( [4.131 ). 

We now assume that \I/^(rsQ,) is analytic, which imphes, as in part (i), that both xiil) 



and XuiQ) ^I's analytic, and hence xI(Q') has a Taylor series 



Xfj, 



(?l,?2,g3) = X]^"A'*^^l'^2)g3: 



(4.47) 



n=0 



where bnniqi, q'2) is analytic. We must now only show that the appropriate coefficients vanish 
in the Taylor series Eq. ( [4.47 ) to produce the Taylor series Eq. ( 2.20 ). Since cr(g) is analytic 
and even in q^, Lemma |l| shows that a{rsa) = cr(5'(i"sa)) is analytic. Therefore the function 
defined by 



^L{qi,q2,u,,u,) = e-^'^(''-''-"-"^)vI/^(gi,g2,«i,«2,e = z) 



(4.48) 



is analytic. Since e = z, we see from Eq. ( [4.15[ ) that the matrix product UVW appearing in 
Eq. (WM) is equal to W^VW. From Eqs. ( ^l2D and ( ^lil) it is evident that W^VWz = z, 
and hence W^VW is a rotation by 6 about the z-axis, where 6 is the rotation angle of V. 
Therefore, 



Pi*^(W^VW)=5„^exp(zm6 



(4.49) 



Thus, Eqs. ( ^4.46[ ) and ([4.48| ) combine to yield 



(4.50) 



where we have dropped the explicit dependence on qi and q2- From this equation and the 
fact that {9, q^) are the usual polar coordinates with respect to the Cartesian coordinates 
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U2), we may apply Theorem |[ In particular, since ^^(mi, ^2) is analytic, Xmi.l'i) has the 
Taylor series given in Eq. (|3.9| ). This implies that the appropriate coefficients in Eq. (f4.47|) 
vanish. 

Next assume that x^(Q')) and hence is analytic and that x^(Q') has the Taylor series 

given by Eq. ( p.20| ). Then from Theorem 0, the function U2) given in Eq. ( |4.5U| ) is ana- 

lytic in Ml and U2 and based on the analyticity of the coefficients a^nl*?!, 92), ^ml""!) ""2, ^i, '?2) 
is analytic in qi and q2 as well. We now rewrite Eq. (^4.46|) as 

< = e"'J2'^t,ime"''xi = e^^'Y.'^t.imK^ (4.51) 

where we have used Eq. ( |4.49| ) in the first equality and Eq. ( [4.50| ) in the second. We have 
already shown that a, U, W, and \E'^ are analytic functions of the mixed coordinates. Thus, 
is an analytic function of the mixed coordinates and hence of the Jacobi coordinates as 
well. QSV. 

V. CONCLUSIONS 

Assuming valid conventions. Theorem |l| solves the problem of determining whether an 
internal wave function for the three-body problem is associated with an analytic external 
wave function. The criteria presented are both necessary and sufficient. We now propose 
several useful extensions of the present work which we plan to consider in future publications. 

First, we observe that Theorem |I] says nothing about the properties of the external wave 
function in the neighborhood of body frame singularities. Our justification for ignoring these 
singularities at present is that they may always be moved away from the region of interest 
by a change of body frame; they may even be moved into the unphysical region of shape 
space. However, in practice it is not always convenient to change body frames or to use a 
body frame which places the singularities in the unphysical region. Therefore, one should 
like to have a set of criteria on the internal wave function which guarantees the analyticity 
of the external wave function in the neighborhood of a body frame singularity. 
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A second desirable extension of the present work would be to develop a set of criteria 
applicable at the three-body collision. Such an analysis would almost certainly involve the 
analysis of body frame singularities mentioned above. This is because the string (or strings) 
of body frame singularities emanates from the three-body collision and no body frame avoids 
singularities at the three-body collision. 

As mentioned in the introduction, the external wave function for three-body problems 
is not necessarily analytic at all points, the collisional configurations of Coulomb problems 
being an obvious example. Yet another useful extension of the present work would be to 
develop criteria on the internal wave function which capture such singular behavior of the 
external wave function. 

Finally, it is natural to try and extend our analysis to four or more bodies. The case 
of four bodies is certainly more challenging than that of three bodies. For the three-body 
problem, we can always pick valid reference orientations ra{q) which are well defined ana- 
lytic functions in the neighborhood of an arbitrary coUinear shape (except the three-body 
collision). For the four-body problem, however, any reference orientation ra{q), a = 1,2,3, 
we choose will not be analytic at any coUinear shape. This is because all body frames for 
the four body problem have singular surfaces emanating from all coUinear shapes [^. Ex- 
tending the analysis of this paper to the four body problem will necessarily require a deeper 
consideration of frame singularities. 
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APPENDIX A: FACTS CONCERNING SO{2) 



We collect some basic facts concerning S0{2) which are necessary for the proof of The- 
orem 0. All irreducible representations (irreps) of SO (2) are one-dimensional and may be 
labeled by an integer m, with — oo < m < oo. The components of a vector in an invariant 
one-dimensional carrier space labeled by m are multiplied by exp{—im6) when rotated by 

e. 

The fundamental representation of 5*0(2), that is the representation by 2 x 2 real or- 
thogonal matrices, decomposes into the direct sum of the m = ±1 irreps, which we denote 
-|-1 © — 1. The two irreducible carrier spaces are spanned by the vectors e± = (1, ±i)"^. The 
d-fold tensor product of the fundamental representation contains (with various multiplici- 
ties) the irreps d, d — 2, d — A, —{d — 4), — (rf — 2), —d. However, exactly one irrep for 
each allowed value of m is totally symmetric. That is, the fully symmetrized tensor product 
of the fundamental representation of 5*0(2) decomposes as 

S 0(+l ® -1) = d ® {d - 2) ® (d ~ A) ® ■ ■ ■ ® -{d - A) ® -{d - 2) ® -d, (Al) 

d 

where ® denotes the tensor product and S denotes the total symmetrization operator. The 
irrep labeled by m is spanned by a totally symmetric rank d tensor tf^ given by 

t^ = 5(e+®---®e+®e_®---®e_), (A2) 

"•^ V ' V ' 

u times V times 

where m = u — v and d = u + v. Explicitly, tj^ transforms as 

ki-.-ka 

where (^m)ji-- id the components of tf^. A final fact we need, which follows readily from 
Eq. (^), is that 

(tt)..... = 1- (A4) 
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